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Abstract:  Adaptive  Hydraulics  (ADH)  is  a  computational  fluid  dynamics 
package  that  solves  the  Navier-Stokes  equations  and  shallow-water 
(depth-averaged  Navier-Stokes)  equations  on  two-  and  three-dimensional 
computational  meshes  for  overland  flow,  sediment  transport,  and  ground- 
water  problems.  The  capability  to  model  the  hydrodynamic  effects  of 
vessels  moving  through  a  two-dimensional  flow  field  has  been  added  to 
ADH.  Using  empirical  relations  developed  by  Maynord  (2000),  the  bed 
shear  stresses  induced  by  a  barge  bow  and  towboat  propeller  may  also  be 
calculated,  which,  in  turn,  can  be  used  to  predict  sediment  transport. 
Guidelines  have  been  developed  for  sufficient  mesh  refinement  near  the 
vessel  and  appropriate  values  for  the  mesh  adaption  parameters  in  ADH. 
To  demonstrate  the  new  modeling  capability,  this  report  describes  in 
detail  a  study  of  the  effects  of  a  vessel  sailing  through  a  stretch  of  the 
Illinois  Waterway  near  Kampsville,  IL.  The  vessel-induced  shear  stresses 
are  also  reported  for  the  Kampsville  study. 
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1  Introduction 

Background 

The  hydrodynamic  effects  of  a  vessel  moving  through  a  flow  field  signifi¬ 
cantly  affect  a  navigation  channel.  The  waves  produced  by  moving  vessels 
increase  the  rate  of  shoreline  erosion  and  affect  other  vessels  moving 
through  the  channel.  Vessel  movement  induces  shear  stresses  on  the 
channel  bed,  which  alters  sediment  transport  and  deposition  patterns. 
Successful  modeling  of  such  effects  provides  a  valuable  tool  for  predicting 
the  impact  vessels  have  on  the  channel  and  other  vessels. 

The  numerical  model  code  developed  in  this  research  followed  the 
methods  that  had  been  used  in  the  original  Upper  Mississippi  River  Study 
(UMRS)  conducted  in  the  mid-to-late  1990s  (Stockstill  and  Berger  1999). 
The  current  study  improved  on  the  former  in  that  the  numerical  scheme 
employed  mesh  adaption  and  incorporated  relations  to  calculate  shear 
stresses  on  the  bed  induced  by  a  moving  vessel.  The  result  is  a  single,  effi¬ 
cient  model  developed  to  do  the  work  that  previously  required  several 
standalone  models.  The  current  research  provides  a  method  of  computing 
the  vessel-generated  bed  shear  stresses  using  the  algorithms  and  data 
presented  by  Maynord  (2000).  The  empirical  shear-stress  relations  given 
by  Maynord  (2000)  rely  on  rules  as  to  the  direction  the  vessel  travels 
relative  to  river  currents.  Rather  than  tuning  the  empirical  parameters  for 
a  particular  channel  configuration  and  vessel  size  and  speed,  a  combina¬ 
tion  of  physics-based  and  empirical  relations  was  deemed  appropriate  to 
solve  the  problem. 

Vessel  effects 

Quantifying  the  effects  of  vessels  sailing  along  a  waterway  is  a  complex 
problem  because  the  flow  is  unsteady  and  the  river  bathymetry  is  nonuni¬ 
form.  These  effects  include:  drawdown  of  the  water  surface,  return  cur¬ 
rents  around  the  vessel,  and  increased  bed  shear  stresses  and  pressures 
that  act  to  re-suspend  the  bed  material  and  perhaps  erode  banks.  Histor¬ 
ically,  studies  to  capture  vessel  effects  required  that  the  flow  first  be 
modeled,  then  the  bed  effects  were  computed,  and  finally,  the  fate  of  the 
bed  material  that  is  re-suspended  was  determined.  A  seamless  modeling 
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system  would  allow  the  hydraulic  engineer  conducting  the  study  to  focus 
more  on  model  results  and  remediation  than  on  the  modeling  process. 

Previous  studies 

Flow  fields  associated  with  moving  vessels  have  been  simulated  using  the 
hydrodynamics  code  HIVEL2D.  This  model  computes  the  vessel  effects  in 
conjunction  with  ambient  conditions  (Stockstill  and  Berger  1999,  2001; 
Stockstill  et  al.  1995).  HIVEL2D  uses  an  unstructured,  static  mesh  wherein 
the  vessel  is  simulated  by  a  pressure  field  that  moves  through  the  water¬ 
way.  The  model  produces  current  and  water-surface  knowledge  such  as 
drawdown  and  transverse  stern  waves.  However,  HIVEL2D  requires  a  fine 
mesh  throughout  the  entire  sailing  path  of  the  vessel.  The  complicated 
meshes  require  long  setup  and  long  computing  times,  and  the  hydrody¬ 
namics  is  not  coupled  with  any  sediment  transport  simulator,  which 
means  the  numerical  model  does  not  directly  include  sedimentation 
effects. 

Current  approach 

The  current  study  addresses  these  issues  using  the  Adaptive  Hydraulics 
(ADH)  flow  solver.  The  shallow-water  module  of  ADH  solves  the  depth- 
averaged  two-dimensional  (2D)  flows  using  serial  or  parallel  processing. 
However,  one  of  the  most  beneficial  attributes  of  ADH  is  the  adaptive 
mesh  technique,  which  gives  it  the  capability  of  refining  or  coarsening  the 
mesh  based  on  error  estimates  during  flow  calculations  (Tate  et  al.  2006). 
Similar  to  HIVEL2D,  the  vessel  effects  are  modeled  in  ADH  by  moving  a 
pressure  field  through  the  flow.  Due  to  mesh  adaption  abilities,  the  code  is 
able  to  improve  upon  the  previous  technique  by  reducing  computational 
times.  Mesh  refinement  can  be  limited  to  the  vessel  vicinity  even  though 
the  position  changes  with  time.  Berger  and  Lee  (2005)  have  documented 
the  mesh  adaption  algorithms  and,  in  particular,  their  application  to 
modeling  moving  vessels.  As  a  vessel  approaches  a  region  within  the 
channel,  the  mesh  in  that  region  is  refined  for  the  calculations  in  the  more 
complex  flow.  As  the  vessel  continues  past  that  point,  the  mesh  is  then 
coarsened.  This  automated  mesh  adaption  reduces  the  cost  of  studying  the 
impact  vessels  have  in  confined  waterways.  Because  the  refinement  indi¬ 
cator  is  determined  by  the  flow  variables,  ADH  will  adapt  any  region  in 
which  the  flow  is  strongly  disturbed.  Generally  this  will  be  in  the  area  near 
the  vessel.  However,  if  the  vessel-induced  waves  are  still  significant  at 
some  distance  away,  ADH  will  recognize  this  and  adapt. 
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Given  input  including  boundary  conditions,  computational  meshes,  and 
vessel  parameters,  ADH  produces  a  resulting  flow  field  and  water-surface 
elevation  by  providing  a  finite-element  solution  to  the  2D  shallow-water 
equations.1  As  the  vessel  moves  along  a  prescribed  sailing  line,  the  effect  of 
its  pressure  field  causes  changes  in  the  flow  field.  ADH  returns  the  com¬ 
puted  shear  stress  as  a  series  of  scalar  values  as  it  does  with  other  output 
variables  (residual,  depth,  and  velocity  components).  The  shear  stresses 
can  be  used  to  compute  sediment  entrainment,  and  the  computed  flow 
field  can  be  used  to  model  sediment  transport. 

Garcia  et  al.  (2002)  and  Maynord  (2000)  have  studied  vessel-generated 
stresses  on  the  bed;  Maynord’s  research  was  selected  as  having  the  most 
appropriate  results  to  incorporate  into  a  numerical  flow  code  to  predict  the 
bed  shear  stresses  that  lead  to  sediment  entrainment  and  thus,  transport. 
The  current  research  is  tasked  with  combining  the  flow  model  that  repre¬ 
sents  the  vessel  effects  with  a  sediment-transport  model  using  a  library 
system.  The  bed  shear  stresses  associated  with  a  vessel  moving  through  a 
waterway  are  quantified  using  empirical  relations  developed  from  lab¬ 
oratory  studies.  The  following  paragraphs  describe  the  essence  of 
Maynord’s  report.  In  various  regions  near  a  vessel,  the  bed  shear  stress  is 
rapidly  increased  due  to  the  presence  of  a  vessel  sailing  over  or  near  the 
particular  region.  These  regions  are  the  focus  of  Maynord’s  study  and  are 
divided  into  areas  where  the  stresses  are  attributed  to  the  vessel’s  bow  and 
to  the  propeller  wash. 

The  approach  Maynord  (2000)  uses  to  calculate  the  vessel-generated  bed 
shear  stresses  requires  knowledge  of  whether  the  vessel  is  upbound 
(moving  in  a  direction  opposite  to  ambient  currents)  or  downbound.  This 
method  was  used  in  the  Upper  Mississippi  River  Study  (Maynord  2000) 
and  by  Garcia  et  al.  (2002)  to  quantity  the  vessel  effects.  Applications 
using  this  method  require  knowledge  about  ambient  current  directions.  A 
2D  flow  model  capable  of  simulating  vessel-induced  bed  shear  stress  and 
thus  sediment  re-suspension  and  bank  erosion  must  rely  on  a  more  gen¬ 
eral  approach  to  quantify  the  shear  stress. 

A  more  practical  method  is  one  that  uses  a  single  empirical  relation  for  the 
bed  shear  stress  in  both  critical  areas,  at  the  bow  and  propeller,  as  func¬ 
tions  of  fluid  velocity  and  vessel  speed.  Employing  this  method  required 


1  Appendix  A  includes  a  listing  and  brief  discussion  of  the  shallow-water  equations. 
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modification  of  the  drag  coefficients  found  by  Maynord  (2000),  but  the 
resulting  algorithm  does  not  require  knowledge  of  vessel  traveling  direc¬ 
tion  (i.e.,  “upbound”  or  “downbound”). 
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2  Modeling  Vessel-Generated  Waves  and 
Currents 

Modeling  method 

Vessel-generated  currents  and  waves  are  modeled  using  the  2D  unsteady 
shallow-water  equations.  These  equations,  presented  in  Appendix  A,  are 
depth-averaged  equations  of  fluid  motion  that  assume  the  pressure  dis¬ 
tribution  is  hydrostatic.  The  shallow-water  equations  are  capable  of 
resolving  the  long-period  vessel  effects  attributed  to  both  shallow-draft 
vessels  (Stockstill  and  Berger  2001)  and  deep-draft  vessels  (Maynord 
2003). 

Investigators  have  calculated  the  effects  of  a  vessel  sailing  through  a 
waterway  using  the  numerical  hydraulic  code,  HIVEL2D.  However,  this 
model  is  computationally  limited  by  longer  setup  and  computational  times 
than  are  now  necessary.  Also,  HIVEL2D  does  not  include  sedimentation 
effects.  The  present  study  moves  the  concept  of  a  moving  pressure  field, 
currently  in  HIVEL2D,  into  the  modern  ADH.  Also,  the  sediment  entrain¬ 
ment  due  to  drawdown  and  vessel  propellers  is  incorporated  into  ADH. 

Previous  studies  have  relied  on  fixed-mesh  methods  to  model  the  vessel. 
The  difficulty  with  this  method  is  that  the  fine  mesh  resolution  needed  to 
define  the  vessel  is  required  along  the  entire  sailing  line.  This  resolution 
requires  more  effort  in  generating  the  finite  element  mesh  and  signifi¬ 
cantly  increases  computational  time.  As  its  name  implies,  ADH  applies 
mesh  adaption  techniques,  which  are  based  on  a  residual  computed  for 
each  time  step,  to  refine  and  coarsen  a  mesh  based  on  flow  complexity. 

Initial  mesh  resolution  and  adaption 

A  major  principle  of  computational  modeling  is  that  a  model  will  converge 
to  the  governing  equations  as  the  resolution  is  increased.  This  level  of 
mesh  detail  is  usually  defined  by  a  mesh  convergence  testing.  With  an 
adaptive  mesh  model,  convergence  testing  is  made  easier  because  mesh 
refinement  is  done  automatically  during  a  simulation. 

Setting  up  a  moving -vessel  model  in  ADH  requires  a  mesh  that  adequately 
describes  the  bathymetry.  Ideally,  the  mesh  should  only  be  fine  enough  to 
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provide  ADH  convergence.  The  initial  solution  accuracy  is  not  important 
at  this  point  since  this  solution  will  be  used  to  determine  appropriate 
adaption  parameters  for  subsequent  runs. 

The  first  step  in  creating  a  mesh  for  vessel  movement  is  to  determine  the 
sailing  line  that  the  vessel  will  travel.  The  sailing  line  may  consist  of  an 
unlimited  number  of  line  segments  and  arcs.  As  the  vessel  moves  from  one 
segment  in  the  sailing  line  to  another,  the  vessel  is  turned  immediately 
about  its  centroid.  The  transition  from  one  segment  to  the  next  should  be 
smooth  to  avoid  any  sudden,  large  turns.  Similarly,  the  radius  of  any  arc 
sections  of  the  sailing  line  should  be  large.  While  the  sailing  line  defines 
where  a  vessel  moves  within  a  flow  field,  the  sailing  line  is  not  required  to 
begin  or  end  within  the  mesh  boundaries.  If  the  sailing  line  does  extend 
beyond  the  mesh  boundaries,  the  vessel  will  sail  “into”  or  “out  of’  the 
mesh  during  a  simulation.  Solution  convergence  is  difficult  as  the  vessel 
traverses  a  mesh  boundary  in  those  cases. 

The  2D  mesh  can  be  described  in  terms  of  longitudinal  and  lateral  element 
sizes.  Here,  “longitudinal”  means  along  the  vessel’s  sailing  line,  and 
“lateral”  means  normal  to  the  vessel’s  sailing  line.  The  primary  region  of 
interest  is  the  area  in  the  immediate  vicinity  of  the  vessel,  the  near  field. 
“Near”  refers  to  a  rectangular  area  centered  on  the  vessel  center  and 
extending  half  a  vessel  width  on  either  side  of  the  vessel  and  half  a  vessel 
width  forward  and  aft  of  the  vessel.  The  vessel  width  is  taken  as  the  width 
of  the  vessel  at  its  widest  point. 

A  test  scheme  was  developed  to  address  the  lateral  mesh  resolution 
required  both  under  the  vessel  and  in  the  near  field  at  the  port  and  star¬ 
board  sides.  The  vessel  was  centered  laterally  on  a  prescribed  sailing  line, 
and  the  number  of  elements  under  the  vessel  was  varied  sequentially  from 
two  to  eight.  Similarly,  the  number  of  elements  in  the  near  region  of  the 
side  of  the  vessel  was  ranged  from  two  to  eight.  Meshing  guidelines  for 
simulating  vessel  movement  in  HIVEL2D  (Stockstill  and  Berger  1999)  fell 
in  within  the  chosen  scheme.  Solutions  were  generated  with  each  test 
mesh.  The  coarsest  mesh  for  which  ADH  converged  indicates  the  level  of 
lateral  refinement  required  for  the  initial  mesh.  ADH  converged  for  each 
test;  therefore,  two  elements  should  span  the  vessel  width  with  one  ele¬ 
ment  spanning  the  near-field  region  on  either  side  of  the  vessel.  This  level 
of  refinement  should  be  used  around  any  possible  sailing  line  of  each 
vessel  being  simulated. 
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The  number  of  elements  that  span  the  vessel  length  is  related  to  the  ele¬ 
ments  spanning  the  boat  width  by  the  element  aspect  ratio  (AR)  defined 
by: 


AR  = 


(1) 


where: 

Le  =  element  length 

We  =  element  width. 

The  element  aspect  ratio  should  be  about  3  or  less.  Previous  studies 
(Stockstill  and  Berger  1999)  have  shown  that  an  element  length  of  one- 
third  of  the  vessel  length  is  sufficient.  These  element  sizes  are  only  sugges¬ 
tions  for  the  initial  mesh  resolution.  If  these  elements  are  too  large,  ADH 
can  refine  them.  Figure  1  shows  the  recommended  resolution  for  the 
coarse  mesh  around  the  sailing  line.  Note  there  are  two  elements  covered 
laterally  by  the  vessel  footprint  and  three  elements  longitudinally. 


One  final  parameter  that  is  somewhat  independent  of  the  actual  physics  is 
the  length  of  the  time  step  chosen  for  the  simulation.  The  time  step  should 
be  set  so  the  vessel  advances  one  element  length  per  time  step.  Thus,  the 
length  of  the  time  step  (At)  is  calculated  as: 
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At 


(2) 


where: 

Vg  =  vessel  velocity  relative  to  ground. 

Figure  2  illustrates  a  vessel  advancing  by  one  element  length  in  one  time 
step.  The  nodes  along  the  sailing  line  are  shown  in  red,  and  the  vessel  is 
moving  from  left  to  right.  Note  that,  at  time  t,  there  are  two  nodes  along 
the  sailing  line  ahead  of  the  vessel  and  one  aft.  At  the  next  time  step,  time 
t+  At,  there  is  only  one  node  along  the  sailing -line  ahead  of  the  vessel  and 
two  aft. 


Note  that,  in  Figure  2,  the  elements  near  the  sailing  line  are  right  triangles. 
Although  the  elements  are  not  required  to  be  right  triangles,  the  water- 
surface  elevation  contours  from  a  simulation  show  a  vessel  footprint  that 
is  much  more  reflective  of  the  actual  shape  of  the  vessel  than  if  the 
elements  are  arbitrarily  shaped. 

Once  an  initial  coarse  mesh  is  created  and  an  initial  simulation  complete, 
adaption  parameters  can  be  set.  ADH  uses  two  user-defined  parameters, 
the  refinement  tolerance  and  the  maximum  number  of  adaption  levels,  to 
control  the  amount  of  refinement  during  a  simulation.  The  refinement  tol¬ 
erance  value  acts  as  a  “trigger”  value  for  ADH  to  determine  the  elements 
that  will  be  refined.  ADH  compares  the  refinement  tolerance  to  the 
residual  calculated  for  each  element.  If  the  residual  is  larger  than  the 
refinement  tolerance,  the  element  will  be  refined;  otherwise,  no  adaption 
occurs.  The  refinement  tolerance  should  be  set  to  1.0  for  the  initial  simu¬ 
lation  because  ADH  scales  the  residual  values  by  the  refinement  tolerance. 
These  scaled  residuals  are  written  to  the  error  output  file  (*_err.dat). 
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Two  different  attempts  were  made  to  define  a  refinement  tolerance  that 
would  only  refine  the  mesh  near  the  vessel.  In  the  first  attempt,  the  largest 
residual  on  the  boundary  of  the  near-field  region  was  used  as  the  refine¬ 
ment  tolerance.  The  largest  residual  was  typically  found  at  the  stern.  This 
refinement  tolerance  value  resulted  in  refinement  within  the  vessel  foot¬ 
print,  but  did  not  produce  enough  refinement  around  the  vessel.  The 
second  attempt  was  to  use  an  average  residual  value  for  nodes  away  from 
the  vessel  to  specifically  not  refine  the  elements  in  the  far  field,  but  to 
refine  the  elements  around  the  vessel.  Since  the  residual  is  drastically 
lower  in  regions  far  from  the  vessel,  a  representative  number  greater  than 
that  residual  should  be  an  adequate  mesh  refinement  tolerance.  Using  a 
far-field  residual  value  as  the  mesh  refinement  tolerance,  ADH  refined  the 
mesh  only  within  the  near-field  region.  Thus,  a  sufficient  refinement 
tolerance  value  was  determined. 

Figure  3  shows  the  section  of  the  mesh  containing  the  vessel.  The  top  pic¬ 
ture  shows  the  initial,  unadapted  mesh,  and  the  bottom  picture  shows  the 
adapted  mesh.  The  adaption  only  extends  one-half  of  the  vessel  width 
laterally  from  the  vessel  whereas  the  entire  vessel  footprint  has  been 
refined  by  at  least  one  level.  The  longitudinal  refinement  is  limited  to 
about  one-third  of  the  vessel  length  ahead  of  the  vessel  and  one-half  of  the 
vessel  length  aft  the  vessel.  Two  levels  of  refinement  are  required  only  at 
the  bow  and  stern.  This  method  of  choosing  a  mesh  refinement  tolerance 
is  recommended  for  vessel  movement  simulations. 

A  mesh  refinement  study  was  conducted  to  determine  the  necessary  num¬ 
ber  of  adaption  levels  required  to  adequately  model  vessel  movement. 
Using  a  refinement  tolerance  chosen  using  the  far-field  node  method, 
analyses  were  conducted  until  the  solutions  converged. 

Plots  of  both  the  transverse  and  longitudinal  velocity  data  for  the  different 
refinement  levels  showed  that  the  solution  converged  with  two  levels  of 
refinement.  Appendix  B  shows  mesh  convergence  plots  for  a  simulation 
performed  to  compare  model  results  to  field  tests,  which  will  be  discussed 
in  more  detail  in  later  paragraphs. 
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Figure  4  illustrates  the  convergence  plots  showing  transverse  flow  veloci¬ 
ties  as  the  vessel  passes  a  location  in  the  flow.  There  is  no  change  in  the 
flow  solution  for  more  than  two  levels  of  adaption,  indicating  that  only  two 
levels  of  adaption  were  required  to  drive  the  residual  to  a  value  less  than 
the  specified  refinement  tolerance.  The  velocity  at  the  other  locations 
showed  a  similar  convergence  pattern,  so  two  levels  of  refinement  should 
be  used  for  a  vessel  simulation. 

In  summary: 

•  The  elements  in  the  vessel’s  near-field  region  should  be  about  one-half 
of  the  vessel-width  wide  and  one-third  of  the  vessel-length  long. 

•  The  mesh  refinement  tolerance  should  be  set  to  an  average  value  of  the 
residual  at  nodes  far  from  the  vessel. 

•  The  number  of  levels  of  refinement  should  be  set  to  two. 
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Figure  4.  Example  flow  velocity  mesh  convergence. 


Vessel  definition 

A  boat  (*.bt)  file  containing  the  vessel  parameters  must  be  defined  for 
vessel  effects  simulations  and  placed  in  the  project  directory.  Appendix  C 
summarizes  the  necessary  lines  for  the  boat  file.  The  data  in  Table  l 
describe  each  line  in  the  example  boat  file  shown  in  Figure  5.  The 
boundary  condition  (*.bc)  file  must  include  a  line  containing  the  “OP  BT” 
card  to  indicate  a  vessel  simulation.  Figure  5  shows  an  example  boat  file 
for  two  vessels  moving  from  opposite  ends  of  a  channel.  (Appendix  D  gives 
an  additional  example  of  a  boat  file.)  The  “vessels”  described  are  three¬ 
wide  by  four-long  barge  tows. 

Initial  vessel  location 

For  simulations  involving  only  one  vessel,  the  vessel  should  be  placed 
completely  within  the  flow  field  at  the  beginning  of  a  simulation  (i.e.,  the 
vessel  should  be  not  placed  so  it  sails  “into”  the  flow  field).  Convergence  is 
much  more  likely  if  the  vessel  is  initially  within  the  flow  field.  The  first 
coordinate  of  the  sailing  line  defines  the  initial  location  of  the  vessel 
centroid,  so  this  coordinate  should  be  far  enough  into  the  flow  field  that 
the  stern  is  within  the  mesh  boundaries.  After  the  vessel  is  placed  on  the 
mesh,  a  simulation  should  be  run  with  a  stationary  vessel  long  enough  to 
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Table  1.  Example  boat  file  card  description. 


Example  Line 

Meaning 

BOAT  2 

Two  vessels  will  be  simulated. 

FDEF11  -119.0  -1.5  2.9 

The  sailing  line  for  Vessel  1  has  one  segment  beginning  at  (-119.0, 
-1.5),  and  the  vessel  begins  with  an  initial  velocity  of  2.9  m/s. 

DR  FT  1  2.74 

The  draft  of  Vessel  1  is  2.74  m. 

BLEN  1257 

The  length  of  Vessel  1  is  257  m.  (Total  length  of  barges  in  tow) 

BWID  132 

The  width  of  Vessel  1  is  32  m.  (Total  width  of  barges  in  tow.) 

PBOW  1  0.1 

For  Vessel  1,  the  ratio  of  the  bow  length  to  the  vessel  length  is  0.1. 

PSTR  1  0.1 

For  Vessel  1,  the  ratio  of  the  stern  length  to  the  vessel  length  is  0.1. 

CBOW  1  0.95 

For  Vessel  1,  the  fraction  of  the  draft  applied  over  PBOW  is  0.95. 

CSTR  1  0.95 

For  Vessel  1,  the  fraction  of  the  draft  applied  over  PSTR  is  0.95. 

PROP  1  2  2.74  6  52  5 

The  propellers  for  Vessel  1  have  2.74  m  diameter  Kort  nozzles 
spaced  6  m  apart.  The  towboat  length  is  52  m,  and  the  propellers 
are  5  m  from  the  stern. 

SDEF  11  0-3080.0-1.5  2.9 

For  Vessel  1,  the  first  segment  of  the  sailing  path  is  a  line  ending  at 
(-3080.0,  -1.5)  m,  and  the  vessel  begins  with  an  initial  velocity  of 

2.9  m/s. 

FDEF2  1-6955.0-1.5  2.9 

The  sailing  line  for  Vessel  2  has  one  segment  beginning  at  (-6955.0, 
-1.5),  and  the  vessel  begins  with  an  initial  velocity  of  2.9  m/s. 

DR  FT  2  2.74 

The  draft  of  Vessel  2  is  2.74  m. 

BLEN  2  257 

The  length  of  Vessel  2  is  257  m. 

BWID  2  32 

The  width  of  Vessel  2  is  32  m. 

PBOW  2  0.1 

For  Vessel  2,  the  ratio  of  the  bow  length  to  the  vessel  length  is  0.1. 

PSTR  2  0.1 

For  Vessel  2,  the  ratio  of  the  stern  length  to  the  vessel  length  is  0.1. 

CBOW  2  0.95 

For  Vessel  2,  the  ratio  of  the  bow  length  to  the  vessel  draft  is  0.95. 

CSTR  2  0.95 

For  Vessel  2,  the  ratio  of  the  stern  length  to  the  vessel  length  is 

0.95. 

PROP  2  1  2.74  6  52  5 

The  propellers  for  Vessel  2  have  2.74  m  diameter  open  wheel 
nozzles  spaced  6  m  apart.  The  towboat  length  is  52  m,  and  the 
propellers  are  5  m  from  the  stern. 

SDEF  2  10  -3744.0  -1.5  2.9 

For  Vessel  2,  the  first  segment  of  the  sailing  path  is  a  line  ending  at 
(-3744.0,  -1.5)  m,  and  the  vessel  begins  with  an  initial  velocity  of 

2.9  m/s. 

ENDDOO 

End  of  file. 
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Boat  2 

FDEF  1  1-119.0-1.5  2.9 

DRFT  1  2.74 

BLEN  237.7 

BWID  1  32 

PBOW  1  0.1 

PSTR  1  0.1 

CBOW  1  0.95 

CSTR  1  0.95 

PROP  1  2  2.74  6  52.5 

SDEF  11  0-3080.0 -1.5  2.9 

FDEF  2  1-6955.0-1.5  2.9 

DRFT  2  2.74 

BLEN  2  237.7 

BWID  2  32 

PBOW  2  0.1 

PSTR  2  0.1 

CBOW  2  0.95 

CSTR  2  0.95 

PROP  2  2  2.74  6  52.5 

SDEF  2  10  -7-3744.0  -1.5  2.9 

ENDDOO 

Figure  5.  Example  boat  file. 


dissipate  the  initial  flow  disturbances.  Depending  on  the  magnitude  of  the 
vessel  velocity,  gradually  increasing  the  velocity  to  its  maximum  may  be 
required. 

Inserting  the  vessel  into  the  flow  field  with  its  full  draft  will  disrupt  the 
flow  field  considerably  and  could  cause  ADH  to  diverge.  If,  for  a  certain 
channel,  ADH  does  not  converge  using  the  full  vessel  draft,  a  vessel  draft 
smaller  than  the  actual  draft  may  be  used  initially.  Using  the  final  depth 
and  velocity  values  in  the  *_dep.dat  and  *_vel.dat  files  as  the  new  initial 
flow  conditions,  the  draft  can  be  increased  and  a  new  simulation  run.  This 
procedure  should  be  followed  until  a  simulation  with  the  full  vessel  draft  is 
completed. 

For  simulations  with  multiple  vessels,  especially  simulations  over  very 
long  periods  of  time,  some  vessels  may  not  begin  a  simulation  in  the  flow 
field  and  will  enter  at  some  time  during  the  simulation. 
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Multiple  vessel  considerations 

ADH  is  capable  of  modeling  multiple  vessels  within  the  flow  field. 

Analyses  with  multiple  vessels  that  move  through  the  mesh  at  different 
times  can  be  performed  without  extending  the  mesh  to  the  appropriate 
starting  location  for  each  vessel.  The  sailing  line  for  each  vessel  can  be 
extended  outside  the  mesh  boundaries.  The  first  segment  in  each  sailing 
line  should  be  set  to  where  vessel  will  enter  the  mesh  at  the  appropriate 
time  while  traveling  at  the  specified  velocity. 

Model  testing 

A  study  was  performed  to  determine  how  well  the  ADH  calculations  com¬ 
pared  to  field  data  of  a  3-wide  by  4-long  barge  train  taken  from  Bhowmik 
et  al.  (1993).  The  flow  velocity  was  calculated  at  five  points  along  a  cross 
section  of  the  river,  corresponding  to  actual  field  gauge  locations.  Table  2 
lists  the  lateral  location  of  each  gauge  from  the  thalweg.  Negative  distances 
are  to  the  left  of  the  sailing  line  looking  downstream  and  positive  distances 
to  the  right. 


Table  2.  Flow  velocity  gauge  locations. 


Gauge 

Lateral  Location  (m) 

071 

-108 

040 

-60 

332 

69 

642 

87 

999 

100 

Figure  6  schematically  illustrates  the  channel  cross-section  showing  the 
gauge  locations  relative  to  the  vessel.  The  model  results  are  depth- 
averaged  values,  so  the  depth  of  the  field  gauges  is  not  taken  into  account. 
Table  3  lists  the  physical  parameters  of  the  model  conditions. 


Figure  6.  Channel  cross-section  showing  gauge  locations  relative  to  the  vessel. 
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Table  3.  Simulation  parameters. 


Total  Vessel  Length  (m) 

237.7 

Total  Vessel  Width  (m) 

32 

Draft  (m) 

2.74 

Vessel  Velocity  (m/s) 

2.9 

Channel  Flow  Rate  (m3/s) 

628 

Tailwater  Elevation  (m) 

4.7 

Manning’s  n 

0.025 

Time  Step  (s) 

20.8 

Levels  of  Adaption 

2 

Figure  7  shows  velocity  histories  comparing  ADH  results  and  field  data. 
Both  the  longitudinal  and  transverse  velocities  are  those  at  Gauge  999. 

The  time  on  this  plot  has  been  translated,  such  that  the  origin  corresponds 
to  the  time  at  which  the  vessel  bow  passes  the  gauge  locations. 


Time  (min) 

Figure  7.  Velocity  comparison  of  ADH  results  and  field  data. 


The  simulation  data  captures  the  low-frequency  wave  behavior  in  both  the 
longitudinal  and  transverse  directions  very  well.  The  error  in  the  timing  of 
peak  values  has  been  attributed  to  the  field  data  not  being  properly  offset 
to  the  time  the  bow  crosses  the  gauges.  Appendix  E  includes  comparison 
plots  of  ADH  and  field  data  at  each  gauge  location. 
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3  Vessel-Generated  Bed  Shear  Stresses 

The  empirical  relations  developed  by  Maynord  (2000)  and  Garcia  et  al. 
(2002)  describe  the  shear  stresses  in  terms  of  cross-sectional  flow  vari¬ 
ables.  General  or  averaged  cross-sectional  values  such  as  flow  velocity  and 
depth  are  not  readily  available  in  a  multi-dimensional  model.  A  particular 
location  in  a  2D  model  has  knowledge  of  local  flow  variables  such  as  the 
flow  depth  and  velocity  only. 

Vessel-induced  bed  shear  stress  relations 

Vessel-induced  bed  shear  stress  relations  were  originally  developed  to 
account  for  a  vessel  traveling  either  upstream  or  downstream.  A  series  of 
equations  has  been  formulated  using  the  Maynord  (2000)  data  for  the 
peak  bow  stress  that  does  not  require  the  up-  or  downstream  condition. 
The  bed  shear  stresses  induced  by  the  vessel  are  broken  down  into  two 
different  parts:  those  induced  by  the  bow  and  those  induced  by  the  pro¬ 
pellers.  Each  shear  stress  relation  follows  the  same  general  formulation 
with  the  shear  stress  directly  related  to  the  fluid  density  and  the  square  of 
a  characteristic  fluid  velocity. 

Bow  effects 

The  bed  shear  stress  induced  by  the  bow  ( Tbow )  has  been  related  to  the 
vessel  velocity  by: 


r 


bow 


(3) 


where: 

p  =  fluid  density 

Vg  =  vessel  velocity  relative  to  the  ground 
Cbow  =  bow  correction  coefficient. 

The  bow  correction  coefficient  (Cbow)  is  an  empirical  function  of  only  the 
depth-to-draft  ratio  as: 
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C 


bow 


0.068 


depth 

draft 


1.29 


Here,  the  depth-to-draft  ratio  is  computed  using  local  values  as: 

depth  _h  +  P 
draft  d 


(4) 


(5) 


where: 

h  =  flow  depth  under  the  vessel 
P  =  pressure  head  induced  by  the  vessel 
d  =  vessel  draft. 

The  pressure  head  (P)  induced  by  the  vessel  is  calculated  as: 

P  =  —  (6) 

P  9 


where: 

p  =  pressure  on  the  water  surface  representing  the  vessel 
g  =  acceleration  due  to  gravity. 

Propeller  effects 

The  propeller-induced  bed  shear  stress  (Tpr0p),  caused  by  the  accelerated 
flow  exiting  the  propellers,  is  calculated  as: 

1 

^  prop  =^Cfsy shear  (7) 


where: 

Cfs  =  skin  friction  coefficient 

Vshear  =  velocity  used  in  the  calculation  of  the  propeller-induced  bed 
shear  stress. 

The  skin  friction  coefficient  (C/s)  of  the  propeller  is  calculated  as: 
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0.01Dprop 


H, 


prop 


(8) 


where: 

Dprop  =  propeller  diameter 
Hprop  =  depth  under  the  propeller. 


The  depth  under  the  propeller  ( Hpr0p )  is  defined  as  the  distance  from  the 
center  of  the  propeller  to  the  channel  bottom  or: 


H 


prop 


(9) 


As  shown  in  Figure  8,  this  relationship  assumes  that  the  propeller  is 
essentially  at  the  water  surface.  This  assumption  is  valid  for  the  barge 
trains  used  to  develop  the  bed  shear  stress  relations  given  in  this  report. 


Vessel 

]  prop 


Figure  8.  Schematic  of  propeller  terms. 

The  velocity  used  in  the  calculation  of  the  propeller-induced  bed  shear 
stresses  is  defined  as: 


▼ 

h 

i  — 

TT 

n prop 

f 

8- 

^ shear 


(10) 


where: 

Vprop  =  maximum  bottom  velocity  from  the  propeller  relative  to  the 
vessel 

Vwake  =  wake  velocity. 
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The  maximum  bottom  flow  velocity  from  the  propeller  relative  to  the 
vessel  is  given  by: 


V  =EV  f- 

v  prop  ^  v  jetJ 


D„ 


prop 


prop 


(11) 


where: 

E  =  empirical  constant  related  to  the  nozzle  type 
Vjet  =  velocity  of  the  jet  produced  by  the  propeller 
/  =  constant  related  to  the  increase  in  jet  velocity  caused  by  the 
propeller  velocity  and  geometry. 

The  constant/is  defined  as: 


where  q  is: 


f  q,  Vjet  *  o  and  q  >  o 
[  o,  otherwise 


q  =  l  —  C 


V. 


Vjet 


H, 


prop 


D, 


prop  ) 


1.5 


(12) 


(13) 


where  values  for  C  are  given  in  Table  4. 


Table  4.  Nozzle  parameters. 


Nozzle  Type 

Do 

E 

C 

P(hp) 

Ai 

A2 

Open  Wheel 

0.7  ID  prop 

0.43 

0.50 

4550 

23.57 

2.3 

Kort 

Dprop 

0.58 

0.25 

3780 

31.82 

5.4 

Note  that  the  factor/cannot  be  negative,  and  therefore  cannot  decrease 
the  maximum  bottom  velocity  from  the  propeller  relative  to  the  vessel.  The 
velocity  of  the  jet  produced  by  the  propeller,  Vjet,  is  given  by: 


V 


jet 


1.13 


D 


(14) 
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where: 

Do  =  the  jet  diameter  at  the  location  of  maximum  contraction  of  the 
jet.  (Table  4  lists  expressions  for  Do.) 

The  thrust  is  solely  dependent  on  the  vessel  velocity  relative  to  the  ambi¬ 
ent  flow  and  the  power  of  the  propeller  engines.  This  relationship  is  char¬ 
acterized  by  Toutant’s  (1982)  equation: 

?;„=a  as) 

where: 


Tprop  =  thrust  per  propeller 
Pves  =  total  vessel  power  (in  hp) 

Vw  =  vessel  velocity  relative  to  the  flow  (in  mi /hr) 

A1}  A2  =  empirical  constants. 

The  vessel  power  and  the  constants  Ax  and  A2  depend  on  the  type  of  nozzle 
on  the  propeller.  Table  4  lists  the  values  for  all  constants  required  to 
calculate  the  propeller-induced  bed  shear  stresses  for  both  open- wheel 
and  Kort  nozzles.  The  maximum  wake  velocity  relative  to  the  ambient  flow 
velocity  is  given  by: 


V  ,  = -0.78V 

wake,  max  g 


depth 

draft 


Each  nozzle  parameter  in  the  shear  stress  model  has  now  been  introduced 
and  defined.  Table  4  lists  the  nozzle  parameter  values  used  in  the  model  as 
reported  by  Maynord  (2000). 

The  wake  velocity  decays  with  distance  from  the  vessel;  however,  the  decay 
rate  for  the  phenomenon  remains  a  research  topic.  As  a  matter  of  com¬ 
pleteness  and  for  simplicity,  a  linear  rate-of-decay  of  wake  velocity  is 
assumed.  The  decay  rate  ( fiwake )  is  defined  as: 
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The  wake  velocity  astern  of  the  vessel  is  given  by: 


=  V. 


wake, max*  wake 


P wake 


(18) 


where: 

Va  =  ambient  velocity  ahead  of  the  vessel. 

Vessel-averaged  pressure  and  velocity 

Vessel-induced  forces  on  the  bed  are  dependent  on  the  depth-to-draft 
ratio.  This  ratio,  defined  in  Equation  5,  requires  a  single  value  for  the 
pressure  head  at  the  vessel.  Multiple  nodes  fall  within  the  vessel  footprint, 
and  each  one  has  a  unique  pressure.  A  single  value  of  pressure  head  is 
needed  to  compute  bed  shear  stresses,  so  they  must  be  combined  in  some 
way.  To  obtain  a  single  value,  the  pressure  head  is  calculated  for  15  dif¬ 
ferent  locations  within  the  vessel  footprint  (Figure  9).  The  pressure  head 
for  the  vessel  is  the  weighted  average  of  the  pressure  head  of  each  of  the 
15  points  and  the  area  of  the  section  of  the  vessel  footprint  in  which  each 
point  lies. 


Figure  9.  Distribution  of  points  used  to  calculate  the  vessel  pressure  head. 
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The  ambient  flow  velocity  is  required  to  calculate  the  vessel  velocity 
relative  to  the  flow,  which,  in  turn,  is  used  in  the  shear  stress  relations 
described  below.  This  ambient  velocity  is  calculated  as  the  flow  velocity  at 
a  point  located  three-fourths  of  a  vessel  length  in  front  of  the  vessel.  The 
velocity  at  this  location,  relative  to  the  vessel,  is  interpolated  each  time 
step.  Care  should  be  taken  during  mesh  construction  to  ensure  that  this 
interpolation  point  never  falls  outside  the  mesh.  That  is,  the  front  of  the 
vessel  should  not  come  within  three-fourths  of  a  vessel  length  from  the 
edge  of  the  mesh.  Figure  to  schematically  shows  the  point  used  to  calcu¬ 
late  the  ambient  flow  velocity  and  its  location  relative  to  the  vessel. 


Vessel  Centroid 

\ 


Ambient  Velocity  Point 


- ► 

Vessel  Direction 


Figure  10.  Location  of  ambient  velocity  calculation  point  relative  to  the  vessel. 


Limitations  of  the  shear  stress  relations 

The  vessel-induced  bed  shear  stresses  routines  are  only  valid  for  barge 
tows  and  towboats,  as  there  are  currently  no  data  to  formulate  vessel¬ 
generated  bed  shear  stress  routines  for  ships.  Also,  the  shear  stress 
relations  must  be  expressed  in  SI  units.  The  stress  relations  described  in 
this  report  give  stresses  in  units  of  dynes/cm2. 


The  trail  representing  the  decay  of  bed  shear  stress  induced  by  the  pro¬ 
peller  is  always  aligned  with  the  vessel.  This  alignment  becomes  an  issue 
with  curved  sections  of  sailing  lines  such  as  when  a  vessel  navigates  a 
bend. 
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Defining  the  propeller  in  ADH 

An  OP  BTS  card1  must  be  added  to  the  boundary  condition  file,  and  a 
PROP  card  added  to  the  boat  definition  file  when  modeling  a  vessel.  Also, 
the  pertinent  propeller  information  must  be  added  to  the  boat  file.  A 
description  of  each  parameter  that  must  be  included  with  the  PROP  card  is 
given  in  Appendix  C.  The  towboat  length  is  also  included  in  the  PROP 
card.  This  length  provides  an  offset  distance  of  the  propeller  induced  shear 
stresses  from  the  vessel. 

Stress  output 

The  stress  values  are  output  to  an  American  Standard  Code  for  Infor¬ 
mation  Interchange  (ASCII)  data  file  that  can  be  read  directly  into  an  XMS 
application  such  as  SMS  or  GMS.2  This  file  can  be  opened  in  SMS  to  show 
shear  stress  contours  just  as  velocity  magnitude  or  depth  contours  are 
generated  from  velocity  (*_ovl.dat)  and  depth  (*_dep.dat)  files.  The 
naming  convention  for  the  near-field  vessel -induced  shear  stress  file  in 
ADH  is  to  append  “_str.dat”  to  the  job  name. 

Bed  shear  stresses  attributed  to  the  flow  field,  including  the  return  cur¬ 
rents,  can  be  computed  using  the  velocity  solution.  These  stresses  can  be 
calculated  in  SMS  (with  the  data  calculator)  using: 

V=fp  CVf  <19> 


where: 

Tf  =  bed  shear  stress  due  to  ambient  and  return  currents 
C  =  coefficient  of  friction  for  the  bed  (taken  as  o.oi  for  a  sand  bed) 
Vf  =  flow  velocity  magnitude  (ambient). 


1  “Card”  refers  to  the  set  of  characters  that  begin  a  line  in  an  ADH  input  file.  In  this  case,  the  line  in  the 
input  file  would  begin  with  “OP  BTS.” 

2  SMS  and  GMS  are  two  software  packages  used  by  the  U.S.  Army  Corps  of  Engineers  to  generate  2D 
(SMS)  and  3D  (GMS)  grids  as  well  as  perform  data  analyses  and  other  post-processing. 
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4  Modeling  Bed  Shear  Stress 

Introduction 

The  vessel-induced  bed  shear  stresses  are  important  because  sediment 
suspension  and  transport  is  determined  by  the  shear  stress  acting  on  the 
bed.  For  a  sediment  particle  to  be  suspended  and  transported  within  the 
flow,  the  shear  stress  acting  on  the  particle  must  be  higher  than  a  critical 
shear  stress  value.  The  shear  stress  induced  on  the  bed  by  a  passing  vessel 
is  a  large  portion  of  the  total  bed  shear  stress,  so  calculating  the  vessel  in¬ 
duced  bed  shear  stress  is  an  important  in  determining  sediment  transport 
patterns. 

Model  setup 

The  bed  shear  stresses  induced  by  the  W.C.  Norman  sailing  on  the  Illinois 
Waterway  at  Kampsville,  IL,  were  calculated.  Chapter  2  describes  the 
details  of  the  vessel  and  the  channel. 

Shear  stress  results 

Figure  11  shows  a  longitudinal  distribution  of  the  bow-,  propeller-,  and 
flow-field-induced  bed  shear  stress  at  an  instance  in  time  during  a  simu¬ 
lation.  These  shear  stresses  are  normalized  by  the  peak  bow-induced  shear 
stress,  and  the  x-axis  shows  the  relative  position  of  the  vessel  with  0.5 
corresponding  to  the  bow  and  -0.5  corresponding  to  the  stern.  The  peak 
bow  shear  stress  occurs  slightly  aft  of  the  bow.  Laboratory  work  by 
Maynord  (2000)  suggests  that  the  peak  bow  stress  occurs  at  a  point  0.04  L 
astern  of  the  bow.  The  propeller-induced  shear  stress  peaks  aft  of  the 
vessel’s  stern  immediately  behind  the  towboat.  Notice  that  the  flow-field- 
induced  shear  stresses,  seen  in  the  regions  greater  than  half  a  vessel  length 
in  front  of  the  vessel  and  two  vessel  lengths  behind  the  barge  train,  are  at 
least  an  order  of  magnitude  smaller  than  those  induced  by  the  vessel. 
Therefore,  the  flow-field-induced  bed  shear  stresses  can  be  ignored  when 
only  the  peak  shear  stress  values  are  of  interest. 
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Figure  11.  Longitudinal  shear  stress  plot  in  the  vessel  vicinity 
along  the  sailing  line. 


Figure  12  shows  the  vessel-induced  bed  shear  stress  contours  at  an  instant 
in  time  as  the  vessel  is  moving  from  right  to  left  through  the  channel.  The 
white  rectangle  represents  the  vessel  outline.  The  contours  show  that  the 
bow-  and  propeller-induced  bed  shear  stresses  are  localized  around  the 
vessel.  The  propeller-induced  stress  is  much  larger  than  that  of  the  bow. 
The  contours  show  that  the  propeller-induced  bed  shear  stresses  decay  but 
remain  aligned  with  the  vessel.  The  stress  contours  are  not  symmetric  with 
respect  to  the  centerline  of  the  vessel  because  the  mesh  is  not  symmetrical. 


Figure  12.  Shear  stress  contours  and  their  location  relative  to  the  vessel. 
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5  Summary,  Conclusions,  and 
Recommendations 

Summary 

Vessel  entrainment  relations  have  been  added  to  the  2D,  shallow-water 
module  of  the  hydrodynamic  finite-element  code,  ADH.  These  relations 
allow  the  user  to  calculate  low-frequency,  vessel-induced  waves  within  a 
flow  field  for  a  vessel  moving  along  a  pre-defined  sailing  path.  The 
required  setup-  and  run-time  has  been  reduced  from  that  required  by 
other  models  because  of  the  adaptive  meshing  capabilities  within  ADH. 
Guidelines  have  been  developed  for  generating  an  initial  mesh  near  a 
vessel’s  sailing  path,  and  suggested  values  for  ADH  mesh  refinement 
parameters  have  been  determined.  A  test  case  was  run  to  compare  the 
calculated  flow  field  to  existing  field  data.  Histories  of  the  calculated 
longitudinal  and  transverse  velocities  compare  well  with  field  data. 

The  capability  to  model  bow-  and  propeller-induced  bed  shear  stresses 
attributed  to  barge  tows  has  been  added  to  ADH.  The  relations  used  to 
describe  the  shear  stresses  acting  on  the  bed  are  empirical  based  on  the 
work  reported  by  Maynord  (2000).  These  relations  have  been  generalized 
so  that  a  single  set  of  equations  defines  both  upbound  and  downbound 
vessel  movement.  Example  plots  of  shear  stress  distributions  near  the 
vessel  are  given.  These  relations  show  the  relative  strength  that  the  bow 
and  the  propeller  have  in  inducing  bed  shear  stresses. 

Conclusions 

Guidance  was  developed  for  setting  up  a  vessel  simulation  using  ADH.  The 
elements  near  the  vessel  should  be  approximately  one-half  of  the  vessel- 
width  wide  and  one-third  of  the  vessel -length  long.  This  ensures  that  two 
elements  fall  laterally  and  three  elements  fall  longitudinally  within  the 
vessel  footprint.  The  mesh  refinement  tolerance  should  be  set  to  a  repre¬ 
sentative  value  of  the  residual  at  nodes  in  the  region  of  the  mesh  far  from 
the  vessel.  Two  levels  of  refinement  are  appropriate.  The  propeller- 
induced  peak  shear  stress  is  several  orders  of  magnitude  larger  than  that 
at  the  bow. 
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Recommendations  for  future  work 

The  shear  stress  relations  are  valid  for  specific  types  of  vessels  and  flow 
conditions.  A  laboratory  study  should  be  directed  at  determining  more 
general  relations  applicable  to  a  wider  variety  of  vessels.  These  studies 
should  include  additional  analyses  of  the  propeller  jet,  which  has  been 
shown  to  induce  the  largest  bed  shear  stress. 
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Appendix  A:  Shallow-Water  Equations 

The  shallow-water  (or  long-wave)  equations  are  a  result  of  the  vertical 
integration  of  the  equations  of  mass  and  momentum  conservation  for 
incompressible  flow  under  a  hydrostatic  pressure  assumption.  The  flow 
depth  (h),  the  depth-averaged  x-component  of  velocity  (u),  and  the  depth- 
averaged  y-component  of  velocity  (u)  define  the  dependent  variables  of  the 
fluid  motion.  If  the  fluid  pressure  at  the  surface  is  taken  as  zero  and  the 
free-surface  stresses  are  neglected,  the  shallow-water  equations  are  given 
as: 


dU  dF  dG 
dt  dx  dy 


(Al) 


where: 


hi 


U  = 


uh 

vh 


(A2) 


F 


uh 

u2h  +  —gh2  —  h*^- 
2  P 

uvh  —  h—^~ 

P 


(A3) 


G 


vh 

uvh  —  h—^~ 

P 

vh2  +—gh2  —h^- 

2  P 


(A4) 
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and 


H 


0 


u^fu2  +  v2 


dy  p  dy 


(A5) 


where: 

p  =  fluid  density 
g  =  acceleration  due  to  gravity 
Zb  =  bed  elevation 
n  =  Manning’s  roughness  coefficient 
P  =  pressure  head. 

The  o’s  represent  the  Reynolds  stresses  due  to  turbulence  acting  on  a  unit 
particle  of  the  fluid.  The  first  subscript  indicates  the  direction  of  the  stress, 
and  the  second  indicates  the  face  of  the  unit  particle  on  which  the  stress 
acts.  Using  the  Bousinesq  approach,  the  Reynolds  stresses  are: 


and 


_  du 

(A6) 

<U*  =  2P  vt  — 

ox 

_  dv 

°»=2pv'% 

(A7) 

du  dv 

(A8) 

°Xy=°yX=PVt  ^  +  ^ 

[dx  dy) 

where: 


vt  =  kinematic  eddy  viscosity. 


ERDC/CHL  TR-08-7 


31 


The  eddy  viscosity  (vf)  can  be  determined  from: 

vt  =  BhV  yff 


where: 

B  =  empirical  constant  that  falls  between  o.i  and  1.0 
V  =  magnitude  of  the  flow  velocity 
/  =  Darcy  friction  factor. 


(A9) 
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Appendix  B:  Mesh  Convergence 


Time  (min) 


Figure  Bl.  Mesh  refinement  results,  longitudinal  component  of  velocity  at  Gauge  040. 


Figure  B2.  Mesh  refinement  results,  transverse  component  of  velocity  at  Gauge  040. 
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Time  (min) 


Figure  B3.  Mesh  refinement  results,  longitudinal  component  of  velocity  at  Gauge  071. 


Figure  B4.  Mesh  refinement  results,  transverse  component  of  velocity  at  Gauge  071. 
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Time  (min) 

Figure  B5.  Mesh  refinement  results,  longitudinal  component  of  velocity  at  Gauge  332. 


Figure  B6.  Mesh  refinement  results,  transverse  component  of  velocity  at  Gauge  332. 
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Time  (min) 


Figure  B7.  Mesh  refinement  results,  longitudinal  component  of  velocity  at  Gauge  642. 


Figure  B8.  Mesh  refinement  results,  transverse  component  of  velocity  at  Gauge  642. 
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Time  (min) 


Figure  B9.  Mesh  refinement  results,  longitudinal  component  of  velocity  at  Gauge  999. 


Figure  BIO.  Mesh  refinement  results,  transverse  component  of  velocity  at  Gauge  999. 
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Appendix  C:  Boat  File  Description 

The  boat  file  card  and  field  descriptions  given  below  in  Tables  C1-C5  are 
current  as  of  July  2007.  These  cards  are  subject  to  future  modifications  as 
the  ADH  vessel  entrainment  model  is  extended.  Current  boat  file  card  and 
field  descriptions  can  be  found  at  http://adh.usace.armv.mil/  under  the  “Manuals” 
(last  bulleted  item):  “Quick  Reference  (Carrillo  and  Hallberg  2005).” 


Table  Cl.  Number  of  vessels. 


Field 

Type 

Value 

Description 

1 

Character 

BOAT 

Card  type 

2 

Integer 

>  0 

Number  of  vessels  to  be  modeled 

Table  C2.  Sailing  line  and  initial  vessel  position. 


Field 

Type 

Value 

Description 

1 

Character 

FDEF 

Card  type 

2 

Integer 

>  0 

Vessel  number  sailing  line  being  defined 

3 

Integer 

>  0 

Number  of  segments  in  sailing  line 

4 

Real  number 

# 

x-coordinate  of  the  vessel’s  initial  position 

5 

Real  number 

# 

y-coordinate  of  the  vessel’s  initial  position 

6 

Real  number 

# 

Initial  velocity  magnitude  of  the  vessel 

Table  C3.  Vessel  draft. 


Field 

Type 

Value 

Description 

1 

Character 

DRFT 

Card  type 

2 

Integer 

>  0 

Vessel  number 

3 

Real  number 

>  0 

Vessel  draft 

Table  C4.  Vessel  length. 


Field 

Type 

Value 

Description 

1 

Character 

BLEN 

Card  type 

2 

Integer 

>  0 

Vessel  number 

3 

Real  number 

>  0 

Vessel  length 
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Table  C5.  Vessel  width. 


Field 

Type 

Value 

Description 

1 

Character 

BWID 

Card  type 

2 

Integer 

>  0 

Vessel  number 

3 

Real  number 

>  0 

Vessel  width 

Table  C6.  Bow-to-length  ratio. 


Field 

Type 

Value 

Description 

1 

Character 

PBOW 

Card  type 

2 

Integer 

>  0 

Vessel  number 

3 

Real  number 

>  0 

Ratio  of  the  bow  length  to  the  vessel  length 

Table  C7.  Stern-to-length  ratio. 


Field 

Type 

Value 

Description 

1 

Character 

PSTR 

Card  type 

2 

Integer 

>  0 

Vessel  number 

3 

Real  number 

>0 

Ratio  of  the  stern  length  to  the  vessel  length 

Table  C8.  Fraction  of  draft  applied  to  PBOW. 


Field 

Type 

Value 

Description 

1 

Character 

CBOW 

Card  type 

2 

Integer 

>  0 

Vessel  number 

3 

Real  number 

>  0 

Bow  draft  ratio 

Table  C9.  Fraction  of  draft  applied  to  PSTR. 


Field 

Type 

Value 

Description 

1 

Character 

CSTR 

Card  type 

2 

Integer 

>  0 

Vessel  number 

3 

Real  number 

>  0 

Stern  draft  ratio 

Previous  studies  have  concluded  that  the  flow  solution  is  not  very  sensitive 
to  the  ratios  described  in  Tables  C6-C9.  A  reasonable  value  for  PBOW  and 
PSTR  is  0.1,  and  0.95  is  appropriate  for  CBOW  and  CSTR. 
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Table  CIO.  Sailing-line  segment  definition. 


Field 

Type 

Value 

Description 

1 

Character 

SDEF 

Card  type 

2 

Integer 

>  0 

Vessel  number 

3 

Integer 

>  0 

Sailing  line  segment  number 

4 

Integer 

0  or  1 

Type  of  segment;  (0  for  line;  1  for  arc) 

5 

Real  number 

# 

x-coordinate  of  the  segment  end 

6 

Real  number 

# 

y-coordinate  of  segment  end 

7 

Real  number 

# 

Vessel  velocity  at  the  segment  end 

8 

Real  number 

# 

If  the  segment  is  an  arc,  the  x-coordinate 
of  the  arc  center 

9 

Real  number 

# 

If  the  segment  is  an  arc,  the  y-coordinate 
of  the  arc  center 

10 

Real  number 

±1.0 

If  the  segment  is  an  arc,  the  direction  of 
turn  (+1.0  for  left  turns  and  -1.0  for  right 
turns) 

The  vessel  velocity  is  specified  at  the  end  of  each  sailing-line  segment.  The 
vessel  velocity  is  varied  linearly  along  each  segment  of  the  sailing  line. 


Table  Cll.  Propeller  parameters. 


Field 

Type 

Value 

Description 

1 

Character 

PROP 

Card  type 

2 

Integer 

>  0 

Vessel  number 

3 

Integer 

lor  2 

Propeller  type;  (1  for  open  wheel;  2  for  Kort  nozzle) 

4 

Real  number 

>0 

Propeller  diameter 

5 

Real  number 

>  0 

Distance  between  centers  of  the  propellers 

6 

Real  number 

>0 

Towboat  length 

7 

Real  number 

>  0 

Distance  from  the  propellers  to  the  stern  of  the  towboat 

Table  C12.  End  of  file. 


Field 

Type 

Value 

Description 

1 

Character 

ENDD  0  0 

Card  type 
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Appendix  D:  Example  with  a  Complex  Sailing 
Line 

The  Kampsville,  IL  study  described  in  Chapter  2  has  a  sailing  line  that  is 
only  a  single  line  segment,  so  the  associated  boat  file  does  not  show  how  to 
setup  a  realistic  vessel  maneuver.  The  sailing  line,  shown  in  Figure  Di,  is 
through  Pool  26  of  the  Mississippi  River.  This  path  is  more  complex, 
having  one  line  segment  and  four  arcs. 


The  red  lines  represent  the  vessel  sailing  line  and  the  black  lines  show  the 
mesh  boundaries.  The  points  labeled  as  letters  in  Figure  Di  are  points  that 
must  be  defined  in  the  boat  file.  For  this  sailing  line,  the  vessel  is  initially 
at  Point  A  and  moves  sequentially  along  the  path  to  Point  F.  Points  G-J 
are  the  centers  of  the  four  arcs  in  the  sailing  line.  The  coordinates  of  these 
points  are  listed  in  Table  Di. 
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Table  Dl.  Sailing  line  coordinates  for  Pool  26. 


Label 

X-coordinate  (m) 

Y-coordinate  (m) 

A 

718713.764 

4313644.615 

B 

716648.176 

4311078.278 

C 

715208.524 

4308511.941 

D 

713424.606 

4305726.525 

E 

708667.494 

4304505.950 

F 

706007.266 

4306258.571 

G 

720024.618 

4307497.203 

H 

708968.097 

4310544.656 

1 

709007.958 

4313059.571 

J 

709016.860 

4307931.471 

The  direction  of  turn  must  be  specified  in  the  boat  file  for  the  four  arcs  in 
the  sailing  line,  Segments  2-5.  Table  D2  lists  the  direction  of  turn  for  each 
arc  in  the  sailing  line  as  well  as  the  corresponding  value  that  must  be 
placed  in  the  boat  file. 


Table  D2.  Direction  of  turn  data  for  Pool  26. 


Segment 

Direction  of  Turn 

Boat  File  Value 

2 

Left 

1.0 

3 

Right 

-1.0 

4 

Right 

-1.0 

5 

Right 

-1.0 

Figure  Dl  shows  the  boat  file  for  this  sailing  line  and  vessel,  a  3-wide  by 
5-long  barge  tow,  drafted  at  2.74  m  and  traveling  at  a  constant  2.4  m/s. 


BOAT  1 

FDEF  1  5  718713.764  4313644.615  2.4 

DRFT  1  2.74 

BLEN  1  297.18 

BWID  132 

PBOW  1  0.1 

PSTR  1  0.1 

CBOW  1  0.95 

CSTR  1  0.95 

SDEF  110  716648.176  4311078.278  2.4 

SDEF  12  1  715208.524  4308511.941  2.4  720024.618  4307497.203  1.0 
SDEF  13  1  713424.606  4305726.525  2.4  708968.097  4310544.656  -1.0 
SDEF  14  1  708667.494  4304505.950  2.4  709007.958  4313059.571  -1.0 
SDEF  15  1  706007.266  4306258.571  2.4  709016.860  4307931.471  -1.0 
ENDDOO 


Figure  D2.  Pool  26  boat  file. 
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Appendix  E:  ADH  and  Field  Data  Comparison 


Figure  El.  Longitudinal  and  transverse  velocity  history,  Gauge  040. 


Figure  E2.  Longitudinal  and  transverse  velocity  history,  Gauge  071. 
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Time  (min) 


Figure  E3.  Longitudinal  and  transverse  velocity  history,  Gauge  332. 


Figure  E4.  Longitudinal  and  transverse  velocity  history,  Gauge  642. 
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Figure  E5.  Longitudinal  and  transverse  velocity  history,  Gauge  999. 
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Appendix  F:  List  of  Variables 


Variable 

Ah  A2 


Ar 


fiwake 


c 


Cb 


ow 


Cfs 

D 


At 


Do 

Dprop 

E 


f 


G 

H 

Hprop 

Le 

P 

P 

Pves 

P 

Tprop 

V 


Tprop 

Va 

Vf 

Vg 

Vjet 


Definition 

empirical  constants 
element  aspect  ratio 

propeller-induced  bed  shear  stress  decay  rate 

empirical  constant 

bow  correction  coefficient 

skin  friction  coefficient 

vessel  draft 

length  of  time  step 

the  jet  diameter  at  the  location  of  maximum  contraction  of  the  jet 
propeller  diameter 

empirical  constant  related  to  the  nozzle  type 

constant  related  to  the  increase  in  jet  velocity  caused  by  the  propeller 
velocity  and  geometry 

Acceleration  due  to  gravity 

flow  depth  under  the  vessel 

flow  depth  under  the  propeller 

element  length 

pressure  head  induced  by  the  vessel 
pressure  underneath  the  vessel 
total  vessel  power  (in  hp) 
fluid  density 
thrust  per  propeller 

bed  shear  stress  induced  by  the  flow-field  velocity 

propeller-induced  bed  shear  stress 

flow  velocity  ahead  of  the  vessel 

flow  velocity  magnitude  (ambient) 

vessel  velocity  relative  to  the  ground 

velocity  of  the  jet  produced  by  the  propeller 
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Variable 


Vprop 

Vshear 


vw 

Vwake 


Vwake ,] 


We 


Definition 

maximum  bottom  velocity  from  the  propeller  relative  to  the  vessel 

velocity  used  in  the  calculation  of  the  propeller-induced  bed  shear 
stress 

vessel  velocity  relative  to  the  ambient  flow  (in  m/s) 
wake  velocity  astern  of  the  vessel 

maximum  wake  velocity  relative  to  the  ambient  flow  velocity 
element  width 
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